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Casimir forces can appear between intrusions placed in different media driven by several fluctua- 
tion mechanisms, either in equilibrium or out of it. Herein, we develop a general formalism to obtain 
such forces from the dynamical equations of the fluctuating medium, the statistical properties of 
the driving noise, and the boundary conditions of the intrusions (which simulate the interaction be- 
tween the intrusions and the medium). As a result, an explicit formula for the Casimir force over the 
intrusions is derived. This formalism contains the thermal Casimir effect as a particular limit and 
generalizes the study of the Casimir effect to such systems through their dynamical equations, with 
no appeal to their Hamiltonian, if any exists. In particular, we study the Casimir force between two 
infinite parallel plates with Dirichlet or Neumann boundary conditions, immersed in several media 
with finite correlation lengths (reaction-diffusion system, liquid crystals, and two coupled fields with 
non-Hermitian evolution equations). The driving Gaussian noises have vanishing or finite spatial 
or temporal correlation lengths; in the first case, equilibrium is reobtained and finite correlations 
produce nonequilibrium dynamics. The results obtained show that, generally, nonequilibrium dy- 
namics leads to Casimir forces, whereas Casimir forces are obtained in equilibrium dynamics if the 
stress tensor is anisotropic. 
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I. INTRODUCTION 

There arc many systems in nature which are sub- 
jected to fluctuations, of thermal or quantum origin. For 
such systems, under certain physical conditions, Casimir 
forces, created by the confinement of fluctuations, ex- 
ist and have been calculated (see, e.g., [TJ). The usual 
way to obtain the Casimir forces uses equilibrium tech- 
niques and is therefore valid only for systems in thermo- 
dynamic equilibrium. This means that the fluctuations 
must satisfy a fluctuation-dissipation theorem that guar- 
antees the existence of an equilibrium state, as discussed 
below. Casimir forces for these systems are calculated in 
the spirit of the original work of H. G. Casimir for the 
electromagnetic case [2J. The method takes as a start- 
ing point the Hamiltonian of the system, from which the 
partition function Z = J exp(— /3H ) is calculated, either 
directly or using functional integration [3J. In the cal- 
culation of the partition function one must take into ac- 
count the boundary conditions, that is, the macroscopic 
bodies which are immersed in the system. The parti- 
tion function of the system will have different values for 
different configurations, e.g., different separations of the 
objects. Once the partition function has been obtained, 
its logarithm provides the free energy F. The final step 
required to obtain the Casimir force is the calculation 
of the pressure as the difference in the free energy when 
the configurations of the macroscopic bodies change (for 
example, changing their position, distance or sizes). For 
instance, in the usual Casimir case of forces between two 
flat parallel plates at separation L, the force per unit area 
is given by F c /A — —dF/dL. 

The second approach also takes as a starting point the 



Hamiltonian of the system. However, in this approach 
the Casimir force is derived not from the free energy but 
from the stress tensor T, which is integrated over the sur- 
face of the macroscopic bodies and then averaged over the 
thermal Boltzmann distribution of the associated Hamil- 
tonian exp(— /3if). The approach based on the stress 
tensor has been taken by several authors [1H7]. In fact, 
both approaches are equivalent and valid for equilibrium 
systems only. The reason is that both are based on prop- 
erties which are only valid in equilibrium situations. The 
former uses the thermodynamic relation for the pressure 
as the derivative of the free energy with respect to the 
volume, and the latter uses the Boltzmann distribution 
function, which is only valid for systems in equilibrium. 

On the other hand, other authors have developed a 
dynamical approach [5j I5U10] . Here the starting point is 
an evolution equation for the considered field(s), sup- 
plemented with a noise source term, so that the evo- 
lution of the field takes the form of a Langevin equa- 
tion. Once this equation is solved, the field is inserted 
into the expression for the pressure and the average over 
the noise is taken. As we will see in the next section, 
if the noise is of internal origin, say thermal or quan- 
tum, this description reduces to the equilibrium one, be- 
cause of the fluctuation-dissipation theorem. However, 
the noise does not necessarily have to be internal but can 
have an external origin [IT . for instance, a system in a 
fluctuating temperature gradient [12], subjected to exter- 
nal energy injection such as vibration |13j or electrically 
driven convection [IT] , light incident on a photosensitive 
medium |15j . or spatially and/or temporally correlated 
noise, as considered in Ref. 5J. Recently, Ref. [16] gen- 
eralized this method to a nonequilibrium temperature 
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gradient. In none of these cases can the equilibrium ap- 
proach be applied. Also, the internal dynamics cannot 
satisfy the condition of detailed balance, and therefore 
the internal noise is not described by the fluctuation- 
dissipation relation. In both cases, it is only possible to 
calculate Casimir forces via the dynamical approach. In 
all these nonequilibrium cases a common feature shared 
with the equilibrium Casimir force is that the origin is 
the limitation of the fluctuation spectrum at large wave- 
lengths. They are, therefore, conceptually different from 
other fluctuation-induced phenomena such as ratchets or 
Brownian motors that act at small length scales. 

The plan of the paper is as follows. We start in 
Sect. II by presenting the Langevin equation subjected 
to a general noise. We stress the differences between the 
cases when the Langevin equation derives from an en- 
ergy functional or not, and discuss the implications of 
the fluctuation-dissipation theorem. Section III derives 
the Casimir force from the stress tensor, while Sect. IV 
calculates the actual Casimir force by substituting the 
solution of the Langevin equation into the stress tensor. 

The subsequent sections V, VI, and VII are devoted to 
the application of the formalism to different physical sys- 
tems and different nonequilibrium conditions, that is, dif- 
ferent ways of violating the fluctuation-dissipation the- 
orem. In particular, Sect. V studies a reaction-diffusion 
system with three types of noise: (1) a noise uncorre- 
lated in space and time, (2) a noise exponentially cor- 
related in time, and (3) a spatially homogeneous noise, 
only fluctuating in time. Section VI is devoted to the 
study of a liquid crystal, with an equilibrium noise, sat- 
isfying the fluctuation-dissipation theorem and therefore 
in an equilibrium situation. We continue with a tempo- 
rally correlated noise and finish the section with a maxi- 
mally correlated noise. Finally, to illustrate the power of 
the method, we apply it to a two-field system where the 
evolution equation is non-Hermitian. Usual approaches, 
based on equilibrium properties, have no applicability in 
this case. We finish with some conclusions. 



II. EQUILIBRIUM AND NONEQUILIBRIUM 
FLUCTUATIONS 

The most widely used tool to study the dynamics of 
fluctuations is the Langevin equations and its related 
Fokker-Planck equation. There is a wide literature on 
this subject, in particular using Langevin equations; see, 
for example, Refs. [T7H2T] . 

Let us consider a linear stochastic differential equation 
for the field <p(r, t), 



d t c/ ) = -Mcf, + ar,t), 



(II.1) 



which is a generalization of the Langevin equation to spa- 
tially extended systems. In this equation, M is an oper- 
ator (usually differential) that can be Hermitian or non- 
Hermitian. The operator does not depend on the field 



notation, we have assumed Langevin equations without 
memory, but the generalization to memory kernels is di- 
rect. The term £(r,i) is a Gaussian noise that represents 
the random or stochastic force acting over the field </>, 
and therefore it is the source of fluctuations for 0. It is 
customary to assume that the noise is Gaussian, and its 
averages are 

(£(M)>=0, (IL2) 
(£(r, t)£(r', t')) = Q6(r - v')S(t - t') - h(r - v')S(t - t'), 

where Q is a Hermitian operator that can contain differ- 
ential and integral terms. Differential terms characterize 
noise of conserved quantities, and integral terms noises 
with spatial correlations. The application of this opera- 
tor to the Dirac delta function produces the spatial cor- 
relation distribution h. The noise here is uncorrelated in 
time, although temporal correlations will also be consid- 
ered below. 



Equation (II.ll admits a solution for an initial condi- 
tion (j>°(r) as 



(r,t) 



-Mtj.0 



4> (r) 



-Mt 



[ dre MT e(r,r). (II.3) 
Jo 



In the limit t —> oo, 0(r, t) reaches a stationary state if 
e~ Mt <fi a —> 0. This implies that the eigenvalues of Ai 
must have positive real parts. 

From the Langevin equation, one can construct a func- 
tional Fokker-Planck equation for the probability distri- 
bution P of the field (p. The technique is standard (see, 
e -g-i [13) j an d its solution (which is not normalizable) is 
a Gaussian of the form 



P[<t>] = 



1 

det/C 



a - / dr<t>Kit>/2 



(II.4) 



where the Hermitian operator K, is the solution of the 
equation [T5] 



MKT 



K,- 1 M + = Q, 



(11.5) 



where Ai + denotes the adjoint of AL The probabil- 
ity distribution P depends both on the matrix A4 and 



also on the intensity of the fluctuations Q via Eq. (II. 5 1 
However, to the best of our knowledge, Eq. (II.5) can 



not be solved analytically, so the operator K. cannot be 
expressed in closed form in terms of M. and Q. 

What happens now if the system is at equilibrium? For 
this case there exists an energy functional F (that can be 
either the entropy [TJJ], a Lyapunov functional [22], the 
Hamiltonian, or a free energy |23j ) which is an integral 
over space of a local functional J 7 , which depends on the 
field <p and its gradients. The evolution equation for <p 
can be obtained by generalization to the continuum of the 
thermodynamics of irreversible process (see, e.g., chap- 
ter VII of [T3] ) • This theory relates the time evolution of 
the fields with its conjugated variables or the so-called 
thermodynamic forces, as 



so the Langevin equation ( II. 1 ) is linear. To simplify 



-£* + £(r,t). 



(II.6) 
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Here C is the (symmetric) Onsager operator, which is 
a generalization to continuum of the Onsager matrix of 
transport coefficients. It is also called the dissipation ma- 
trix. The second law of thermodynamics requires that 
the real parts of the Onsager operator eigenvalues are 
positive in order to guarantee local increase of entropy. 
The fields $ appearing in (II. 6) are the conjugated vari- 
ables of <b and can be derived from the functional F as 



SF 



(11.7) 



If the fluctuations are small and the system is far from 
a phase transition, we can assume that SF/S(f> is linear in 
the field </>. This implies that F\<p\ is bilinear in the field 



F[4>] 



(II.; 



and the equation above defines the operator Q. Although 
the definition of J- may not be unique, Q is unique [23] , 
and Q must be positive definitive in order to guarantee 
the existence of a minimum of the free energy. It can 
also be chosen to be Hermitian, because the antisym- 
metric part does not contribute to the free energy. In 
this way, the corresponding Langevin equation is linear 
in d>. Therefore, we can write that 



SF 



(II.9) 



the evolution equation for 



Combination of Eqs. (|II.6|) and (II. 9 1 allows us to write 

as 



8t4> = -cg<i> + s(T,t). 



(11.10) 



If the system is at equilibrium, the well-known 
fluctuation-dissipation theorem |25| imposes that the in- 
tensity of the noise, given by <2, must be related to the 
Onsager operator C by 



Q = k B T(C + C^ 



(11.11) 



Equation (11.10 1 is formally equal to Eq. (II. 1), with the 
operator M given by M = CQ. As both Q and (C + C + ) 
are Hermitian and definitive positive, it can be shown 
that the eigenvalues of M. have positive real parts, even 
though M. can be non-Hermitian or undefined (as in 
the case of the linear hydrodynamic equations) (see [TH] , 
chapter V) [26] . In both equilibrium and nonequilibrium 
dynamics we will assume that the real part of the spec- 
trum of M is strictly positive, i.e., there are no neutral 
modes as happens when there is continuous symmetry 
breaking [57] or critical phenomena [28 , [2S] . 

In equilibrium the fluctuation-dissipation relation has 
drastic consequences for the solution of the Langevin 
equation associated to Eq. (11.10). The equation (II. 5) is 
now written 



k B T(£ + £+), 



(11.12) 



which admits the solution K, 
into Eq 



j3Q. Once substituted 
the probability distribution is given by 
the exponential of the functional F multiplied by /3 — 
(ksT)^ 1 . More precisely, 



(II.4I 



PI 



1 



-PF[4>] 



(11.13) 



where Z is the partition function or the normalization 
constant of P. Given this probability distribution P[<j>], 
we can now calculate the average of any dynamical vari- 
able A((f>) as 



(A) = / dMWPb 



(11.14) 



In particular, the average of the functional F can be cal- 
culated as 



(F) 



dlnZ 

d/3 ' 



(11.15) 



Equations (11.13) and (11.15) are only valid for equilib- 



rium systems, for which an energy functional exists and 
the fluctuation-dissipation theorem is valid. However, if 
the system is out of equilibrium, the probability distribu- 
tion is not the exponential of F and therefore its average 
is not given in terms of the partition function Z. 



III. CASIMIR FORCES FROM THE AVERAGE 
STRESS TENSOR 

How is this discussion related to the calculation of 
Casimir forces? The Casimir force is normally calculated 
for equilibrium situations, that is, when the noise is of 
thermal origin and the fluctuation-dissipation theorem 
is satisfied. One way to calculate the Casimir force is by 
evaluation of the stress tensor T. From the functional F, 
the stress tensor is calculated as [37J 



I; 



T = 



,F - V, 



2V; 



(III.1) 



which allows the definition of the symmetric bilinear 
stress tensor operator T ■ For isotropic systems, the local 
stress is simply given by the diagonal components of the 
stress tensor, or by one-third of its trace. Usual forms 
of T are A(^(r) 2 times the identity matrix or a tensorial 
product of gradient as in liquid crystals, but it can be 
also nonlocal, as in 30J . 

Because of the intrinsically fluctuating nature of the 
fields, the stress tensor has to be averaged over th e ran - 
dom fields £(r, t) or the probability distribution (II.4). 
Once we have the averaged stress tensor, the Casimir 
force over a body of surface S is obtained as 



(T(r)) • hdS, 



(III.2) 
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where the integral extends over the surface of the embed- 
ded bodies and the vector n is a unit vector normal to 
the surface, pointing inward the body. 

As in the original Casimir calculation, the geometry 
that will be considered throughout this paper consists 
of two parallel, infinite plates, perpendicular to the x- 
axis, separated by distance L x (Fig. [I]). In this geometry, 
the Casimir force per unit area on the plates is then the 
difference between the normal stress on the interior and 
exterior side, where the latter is obtained by taking the 
limit L' x — > oo. The force per unit area on the left plate 
is 



F c /A 



(T xx (x = Q;L x ))- lim (T xx (x = 0; L' x )) 

L'—>O0 



(III.3) 

The interpretation is that, if Fc/A is negative, the plates 
repel each other, while if it is positive, an attraction be- 
tween the plates appears. 



o 



x 



FIG. 1. Parallel-plate geometry used to compute the Casimir 
force. The system is confined between plates located at x = 
and x — L x . Additional plates are located at distances L' x 
from these plates, and finally the limit L' x — > oo is taken to 
mimic an infinite system. 



Let us discuss Eq. (III. 2) for equilibrium and noncqui- 
librium situations. In the former case, i.e., a system in 
equilibrium, the average of the stress tensor can be taken 
in two ways: as an average over the probability distribu- 
tion given by Eq. (11.13), or as an average over the fluc- 
tuating term £(r, t). Equilibrium thermodynamics guar- 
antees that both averages are the same. In contrast, in a 
system out of equilibrium, we are left with one option, the 
average over the noise £(r, i), because JC cannot be ob- 
tained in general. As mentioned, the system can be out of 
equilibrium if the fluctuation-dissipation relation is not 
satisfied. In this case, there still exists a functional F 
(from which the Langevin equation is constructed), and 



the stress tensor can be defined via Eq. ( III. 1 ). Then, the 
average in (III. 2 1 has to be taken over the noise. Finally, 



a more complex situation is when the Langevin equation 
is in its most general form, i.e., Eq. (11.1), without M. 



deriving from an Onsager matrix and a functional F. In 
this case, the stress tensor cannot be constructed from 



Eq. ( III. 1 ) , and one must appeal to other considerations 



in order to construct a stress tensor. One can use a micro- 
scopic analysis of momentum transfer, kinetic theory, or 
invoke, for instance, the existence of a hydrostatic pres- 
sure from which the Casimir force can be derived. We 
will assume, therefore, that it will always be possible to 
build the stress tensor operator IT- 



IV. COMPUTATION OF CASIMIR FORCES 

In this section we will develop a formalism, valid for 
both equilibrium and nonequilibrium systems, that al- 
lows us to compute the average stress tensor and there- 
fore the Casimir force. We will assume that the dynamics 
close to the stationary state is described by the dynam- 



ical equation (II.l), where the noise term is assumed to 



be Gaussian with vanishing mean. We assume that the 
noise has temporal and spatial correlations, but no cross- 
correlations, 



(Z(r,tW,t'))=h(r-r')c(t-t'). 



(IV.l) 



Note that we have assumed a dynamical model whose 
deterministic part is local in time (no memory) but that 
the noise can have some memory. This possibility is 
not allowed by the fluctuation-dissipation theorem, and 
therefore the system is automatically put out of equilib- 
rium. A necessary condition to recover equilibrium is a 
local correlation in time, although this condition is not 
sufficient, as shown in Sect. [IT] 

To solve (II.l) we construct the left and right eigen- 



value problems of M with the appropriate boundary con- 
ditions over the immersed bodies. Although we will con- 
sider the case of two parallel plates, the formalism devel- 
oped in this section is completely general. The left and 
right eigenvalue problems read 



Mf n (r) 
M+g n (r) 



Mn/n (r) , 
V*n9n (r) , 



(IV.2) 
(IV.3) 



with the boundary conditions provided by Ai (which are 
the same as those of C if the dynamics derives from a 
free energy functional). The left and right eigenfunctions 
are orthogonal under the scalar product, i.e., (g\f) = 
J dr<7*(r)/(r); that is, under appropriate normalization, 
(<?n|/m) = 5 nm . We can project the field and the noise 
over the left eigenvalues 

«/»(r,t)=^0 n (t)/„(r), £(r,t) = 5^ n (t)/ n (r), 

n n 

(IV.4) 

where cj) n (t) = (g n \(j)(t) ) and f n (i) = (g n \£(t)). By insert- 
i ng t hese expressions ( |IV.4 ) into the evolution equation 
(II.l) we get the evolution equation of each mode <f> n (i) 
as 

OtMt) = -Mn0n(t) + £»(*). 0V.5) 

This equation can be solved analytically as 



0n(t) 



-fJ'n.t 



(IV.6) 



n (O) + / e"» T £ n (r)dT 
Jo 

The first term e~ f4 ™ 4> n (0) is a transient term that van- 
ishes for times longer than t 3> Rc ^ ^ , so that the aver- 
age of each mode over the noise £ is zero in this limit. 
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To compute the average stress tensor at each point, we 
need to compute (T]</>, 4>, r]). Expanding on the eigen- 
value basis and using that = <j>* we get that 

(T(r,t)) = J] (^n(i)CW) T nm (r), (IV.7) 



where T„ m (r) = T[/„, /*„,r]. 

The cross-average of the mode amplitudes is ob- 
tained from (IV.6) and in the stationary regime [t 3> 
l/Re(/x„), l/Re(/i m )] can be written as 



that /i nm = k B T(fi n + /4J (fif„|0 1 .9m), and therefore 
the equilibrium average stress tensor simplifies to 



Jo Jo 



(fn(n)C(^)> • (IV.8) 



Therefore, we need to calculate the correlation of the n 
and m components of the noise 

(Cn(n)C( r 2)) = /rfri l dY 2 g* n {r 1 )g m {r 2 ) {i{Y ll T 1 )^{Y2,T 2 )) 



, , , , (IV.9) 

Substituting Eq. pVT} into pV\9| and pV8| , it is found 
that 

lim (4> n (tW m (t)) = K m C ^ n) + C{ f m \ (IV.10) 
where 

ftnm = J dr 29n( r l) h ( r l~ r 2)g m (r2) = (g n \Qg m ) 

(IV.ll) 

and c is the Laplace transform of c. 

Finally, the local average of the stress tensor in the 
stationary regime, where transients have been eliminated 
and the value of the stress tensor is independent of time, 
is given by 



(T(r)) 



g (Mn) + c(/4J h 



„mTnm(r). (IV. 12) 



This expression is generally divergent when summed over 
all eigenfunctions. This divergence comes from the high- 
est eigenvalues (corresponding to small wavelengths) and 
is due to consider the mesoscopic dynamics given by 
Eq. (II.ll, valid for all wavelengths. However, it is only 



valid above a certain minimal distance (the atomic or 
molecular length, for example). There are some tech- 
niques to avoid this divergence. For instance, a short- 
wavelength cutoff could be introduced as in Ref. [52"] . 
but here we will use regularization techniques similar to 
the Riemann zeta function used in the electrodynamic 
case [2]. 

Using the previous expression, the conditions under 
which Casimir forces exist in an equilibrium system can 
be deduced. As mentioned above, if the dynamics is local 
in time, the fluctuation-dissipation theorem implies that 
the noise terms must not have memory either, therefore 
c(/i) = 1/2. Also, the equilibrium relation (11.11 ) implies 



.{v)) = k B TY J {9n\Q- 1 g m )T nm {Y). (IV.13) 



If the free energy functional depends only on cf> but not on 
its derivatives, the stress tensor operator turns ou t to b e 
isotropic and is given by T„ m (r) = f n Qf^J 3X 3 flllll| ). 
Then, thanks to the completeness of the basis, the stress 
tensor can be further simplified to (T oq (r)) = ksT5(r). 
This expression, once properly regularized, gives a stress 
that is independent of system size; that is, the stress is 
not renormalized by the fluctuations in a size-dependent 
way and therefore no Casimir force can be developed. 
On the contrary, if the stress tensor is not isotropic, as 
in the case of liquid crystals, the result is not trivial and 
Casimir forces can develop, as shown in jSTj . 

All these equations provide expressions for the average 
fields and fluctuations, expressed in terms of the eigen- 
values and eigenvectors of the problem, which encode the 
information of the evolution equation together with the 
boundary conditions. 

To summarize, in this section we have proven that the 
Casimir force over a body is given by 



'c 



E 



T nm {r) ■ hdS, 



M" + Mm JS 

(IV. 14) 

which is the main result of this paper. It shows how to 
derive the Casimir force from the dynamical equations for 
the field <j) subjected to any kind of noise. It is obtained 
by diagonalizing the evolution operator of the field, and 
projecting the noise correlation and the stress tensor over 
the set of eigenfunctions. This approach provides the 
Casimir force for both equilibrium and nonequilibrium 
systems. 



Equation (IV. 14 1 shows the well-known nonadditive 
character of the Casimir force: neither the eigenvalues 
nor the eigenfunctions for different boundary conditions 
are easily related. They cannot be written as a sum of the 
eigenvalues and eigenfunctions of each different problem. 

The rest of the paper deals with applications of 
Eq. (IV. 14) to different physical systems, in both equi- 



librium and nonequilibrium situations. 



V. REACTION-DIFFUSION SYSTEMS 

To show how this formalism works, we calculate the 
Casimir pressure between two plane, infinite plates sep- 
arated by distance L x immersed in a medium described 
by a quadratic free energy 



F[4>] 



dr/ o 2 (r)/2. 



(VI) 



The multiplicative constant /o can be absorbed into </>, 
and we will eliminate it in what follows. The dynamics 



G 



is described by two transport processes: relaxation and 
diffusion; that is, the Onsager operator is C = X — DV 2 , 
where A and D are the transport coefficients (and con- 
sequently, positive) associated with the two irreversible 
processes of relaxation and diffusion, respectively. The 
resulting equation is 



90 
dt 



-\(j>+ DV 2 (t> + £,(r,t). 



(V.2) 



Fluctuation-dissipation is satisfied if the noise is delta 
correlated in time and the space correlation function is 



necessary and other procedures are available. There are 
various regularization methods that allow the isolation 
of the divergent term that is independent of the plate 
separation and therefore cancels out in the computation 
of the Casimir force. The regularization method used in 
this manuscript is based on the Elizaldc function detailed 
in the Appendix. 

To obtain quantitative predictions, we consider specific 
cases for the noise correlations. 



h{v) = 2k B T{X- DV 2 )5{v). 



(V.3) 



A. Uncorrelated noise in time and space 



As the energy functional for this system is cj> 2 /2, with- 
out terms with spatial derivatives, the stress tensor is 
identical to the local energy functional. Also, the dy- 
namic operator is Hermitian, implying that eigenvalues 
are real and that there is no need to distinguish between 
left and right eigenfunctions. In order to obtain Casimir 
forces the appropriate boundary conditions are of Neu- 
mann type, as Dirichlet boundary condition would imply 
trivial vanishing forces. 

We need to solve the eigenfunction problem for the 
spatial part of the dynamics in order to calculate the av- 
erage of the fields that will lead to the Casimir pressure 
over the plates. So, we hav e to solve the eigenfunction 
problem given by Eq. (IV. 2) with A4 = X — DV 2 obeying 



Neumann boundary conditions (no-flux boundary condi- 
tions), d x (/)(0,y,z) = d x (j)(L x ,y,z) = 0. 

The normalized eigenfunctions are characterized by 
three indices n x , n y , and n z , denoted as a whole by n, 
and their form is 



/n(r) 



1 e ik|| 



if n x — 

/„(r) = y/± cos(k x x)e^« if n x > 1. (V.4) 

Here, = yy + zz and ky = k y y + k z z. The eigenvalues 
are 



D ki 



k 2 ) 



(V.5) 



— n k 



j^-riy, and k z 



2tt , 



where k x 

0, 1, 2, . . . and (n y ,n z ) eZ 2 . The quantity fc^ 1 = y/D/X 
is the characteristic correlation length of the system. 
The average stress is then 



(l«(r)) 



2D ^ 



1 L nm 



2k 2 



/n(r)/*(r). 



. (V - 6) 

This expression needs to be regularized, otherwise it is 
divergent. The divergence, as explained in 32 ], is due to 
the application of the mesoscopic model (V.2) up to very 



large wavevectors. Conceptually, the stress could be reg- 
ularized by considering generalized hydrodynamic mod- 
els valid for high wavevectors leading to finite stresses, 
but as the Casimir forces have their origin in the limita- 
tion of the fluctuation at small wavevectors, this is not 



We first consider the case of a noise with vanishing 
correlation time and length, and intensity T, i.e., 



(^r,^(r',t0) = r5(r-r')«5(t-O- 



(V.7) 



This noise correlation, without the — V 2 (5(r) term, au- 
tomatically puts the system out of equilibrium, as Q ^ 
C + C + . The addition of such a term would have led 
to a stress that was independent of plate separation, not 
producing a Casimir force [9]. Therefore, we consider 



the effect of the nonequilibrium noise (V.7) on Casimir 
forces. In this case, h nm (c(^, n ) + c(/i m )) = TS nm and the 



double sum in ( V.6 ) is reduced. On the surface of a plate 
and applying the limit L y ,L z — > oo, the stress is given 
by 



CMO)) 



r 



\q^l x d Jk 



E 



p poo 



k 2 + k 2 + k 2 



8nL r D 



/ dkk g 



(V.8) 

where polar coordinates in the y- and z-components have 
been used. Note that the original sum over n x in (V.6) 



runs only over N, but the form of the normalizations of 
the eigenfunctions (V.4) allows extension of the sum over 



Z with a prefactor of 1/2. 

This expression is divergent, so it must be regularized 
In order to do so, we use the Chowla-Selberg expres- 
sion shown in Eq. (A.3). The parameters are s = 1 
a = tt/L x , and u 2 = k z 



fc 2 . The first term in the sum 



of (A.3) equals L x j ^Jk 2 + k 2 , which combined with the 



prefactor in Eq. ( V.8 ) yields a term which is independent 
of L r , and therefore its contribution to the stress tensor 



cancels in virtue of Eq. (III.3). It must be remarked that 



the size-independent term is actually divergent if the con- 
tinuous model is assumed to be valid for any wavevector. 
Then, we are left with the infinite sum of modified Bessel 
functions K 1 / 2 - This sum can be performed analytically, 
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with the result 



F c /A 



r 

4wD 



dk- 



k 1 

v/fc 2 + k 2 

rfc ln(l - e - 2koL °>) 



- 1 



8ttD 



(V.9) 



Let us note that, because the divergence was eliminated, 
we could have interchanged the integral with the summa- 
tion of the modified Bessel functions to obtain the same 



result. Equation (V.9) shows that the Casimir force di- 



verges if the correlation length tends to infinity, i.e., if 
k — » 0. This result was obtained in [35] using a regular- 
izing kernel technique. 



B. Temporally correlated noise 

We next consider the case of a noise that is delta cor- 
related in space but has exponential correlation in time 

<£ (r, t) £ (r', 0) = TS (r - r') f 1 + ^ e --l*~*'l , (V .10) 



where the factor (1 + |) allows both the white noise limit 
(a — > oo) and the quenched noise limit (a — > 0) to be 
taken. Again, the delta correlation in space leads to a 
term S nm that eliminates one summation in the stress at 
the plate, which is then given by 



(T ra (0)) 



1 



1 



If a > 0, we can factorize the quotient as 



(V.ll) 



Taking the limit L y ,L z — ¥ oo and using polar coordi- 
nates, 



(TWO)) = 



AirL r D 2 



(V.15) 

\ -2 



k -\- kr\ 



Although this expression is finite and does not require 
a regularization procedure, the size-independent contri- 
bution can be eliminated using the same regularization 
procedure as before, using Eq. (A.3) with s — 2. The 
result is 



F c /A 



4irD 2 



dk- 



D-i 



(V.16) 



where ui = \Jk 2 
obtain 



After carrying out the integral, we 



F c /A = 



4irD 2 k e 2k ° L * 



(V.17) 



We remark that this system is not dynamically fluctu- 
ating, because the noise is quenched and the transients 
have been eliminated. Nevertheless, it creates a Casimir 
force whose origin is the same as previously considered 
in the sense that the presence of the second plate limits 
the spectrum of possible fluctuations, and therefore the 
renormalized stresses on the two sides of the plate are 
different. 



C. Maximally spatially correlated noise 



(T ra (0)) 



1 



2aV 



E 



i 



i 



(J-n [In + a 



(V.12) 



with 



— k 2 + k 2 + 

stress is the difference between the Casimir stress of two 



k 2 as before. 



We note that this 



systems with a white temporal noise (V.8) of intensity 



(1 + a)r/o, the first one with fcp = -A and the second one 



with k\ = 



(T xx (0)) 



Then, the stress on the plate is 



r i 



4:airDL x 



In 



-2k L x 



1 



(V.13) 



The Casimir force per unit surface on the plate is given 
just by this expression, because the stress on the un- 
bounded side vanishes [as shown by taking the limit 
L x — >• oo in Eq. (V.13l]. Finally, we can reobtain the 
white noise limit if a — > oo. 

The case a — > corresponds to the quenched limit, 
where static sources of noise are randomly distributed in 
space. The average normal stress at the wall is 



(V.14) 



As a final case, we consider the situation in which the 
medium is perturbed externally by a spatially homoge- 
neous noise, with vanishing correlation time. This could 
be the case when a rapidly fluctuating external field is 
applied to the medium. 



(V.18) 



Applying the same computation procedure as in the other 
cases, the average local stress on each side of the plates 
is 



(TWO)) = 



r l 

2L>2fc2 



r 

4A' 



(V.19) 



which is independent of the plate separation. Therefore, 
the Casimir force vanishes in this case. 



VI. LIQUID CRYSTALS 

The existence of Casimir forces in liquid crystals has 
been known for some time now [31] . In this section we 
apply the formalism presented in Sect. ITT] to a nematic 



crystal, obtaining the known Casimir force for an equi- 
librium situation, and expressions for the force for some 
nonequilibrium conditions. The free energy functional of 
a nematic liquid crystal |33| can be written in terms of a 
planar field <fi as 



F = dr 



2 ' 



H-2 



(VI.1) 



where we have assumed that the director vector is writ- 
ten in terms of the field 4> as n = (sin cf>, 0, cos 4>) , together 
with the one-constant approximation (proportional to 
K2). The first term in Eq. (VI.l) comes from a mag- 



netic field directed along the z-axis, whose intensity is 
absorbed into K\. 

The simplest dynamical model is obtained with a sin- 
gle relaxational transport coefficient, with the Onsager 
operator C = A, leading to 



d(b 

dt 



(VI.2) 



which is identical in form to (V.2|, but with three main 
differences: the form of the fluctuation-dissipation rela- 
tion to be in equilibrium, the stress tensor, and the pos- 
sible boundary conditions that produce Casimir forces. 
Fluctuation-dissipation is realized, according to Sect.pT} 
if the noise satisfies 

(C(r, *)£(r', t')) = 2k B TXS(r - r')S(t - t'), (VI.3) 

i.e., is purely nonconservative. 

Due to the presence of the gradient terms in the free 
energy functional (VI.l I, the stress tensor is not isotropic, 



and therefore even in equilibrium Casimir forces can ap- 
pear. Using Eq. ( III. 1 ) the xx component of the stress 
tensor is 



11 xx — 



'21 

2 



dy 



K2 

2 \dz 



K2 

2 {dx 

(VIA) 

It is then possible to develop Casimir forces by im- 
posing either Dirichlet or Neumann boundary condi- 
tions. Dirichlet boundary conditions are equivalent to the 
strong anchoring conditions, i.e., = over the surfaces, 
and will be the case studied here. Casimir forces with 
Neumann boundary conditions can be easily extracted 
from the Dirichlet ones. 

In this case of Dirichlet boundary conditions, the eigen- 
functions of the operator M. = A [«i — K2V 2 ] are given 

by /n(r) = y y sin (k x x) e lk H' r " , with eigenvalues 



\K 2 {k 2 x 



where k x — -f-n x 



k y — 



2tt , 



K = 



2tt 



11 z 



(VI.5) 

and ko = 



with indices n x GN and (n y ,n z ) G Z 2 



Because of the boundary conditions, the xx component 
of the stress tensor at the plates is simply given by 



T xx (0) 



K2 

2 



c):v 



(VI.6) 



As in the case of the reaction-diffusion system, we will 
consider different types of noise correlations that, as will 
be shown below, produce Casimir forces of different char- 
acter. 



A. Uncorrelated noise in time and space 

We consider an uncorrelated noise as described by 



( V.7 1. This case can be considered as in equilibrium with 



a temperature given by T = 2fcgTA. Again the double 
sum in Eq. (IV. 14 1 can be reduced, and the stress tensor 



on the surface of a plate is given by 



(T XI (0)) 



2VA 

r> 

Applying the limit L yi L z 



E 



fc2 + kl 



(VI.7) 



00, we obtain 



dk, 



E 



n^zjW) +fc2 + fc2 + Ag 



(VI.8) 



Using polar coordinates and regularizing the resulting 
expression using Eqs. (A. 6) and (A.3| with s = 1, 



F c /A 



r 

'4^A 



2 + k 2 o 



dkk 



16ttAL3 
+2^L 2 Li 1 (e- 2fcoi 



(VI.9) 

Li3 ( e -2fc L x ^ + 2k L x U 2 ( e - 2kaL *) 



(VI.10) 



where Li s (z) = X)^=i ^ s polylogarithm function. 
At distances long compared with the correlation lenght, 
that is L x ^> Icq 1 , the force decays as 



F c /A 



Vk 2 

1 K _ e -2k L x 



8ir\L 7 



(VI. 11) 



In the opposite limit, when the plates are at a distance 
much smaller than the correlation lenght, or L x <C , 
the force is 



F c /A 



r C(3) 

16ttA LI ' 



(VI.12) 



This result has already been obtained in the context of 
liquid crystals in [5] if we use the fluctuation-dissipation 
theorem. It is also the high-temperature limit of the 
electromagnetic Casimir force between two perfect metal 
plates i34j (described by Dirichlet boundary conditions) , 
and by using the classical limit of the fluctuation- 
dissipation theorem. 
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B. Temporally correlated noise 



C. Maximally spatially correlated noise 



We consider the temporally correlated noise described 
in Eq. (V.10). By using the eigenfunctions of the Dirich- 



let problem, the stress tensor over a plate takes the value 

« 2 r(i + f)v- i kl 



(Txx(0)> 



E 



(VI.13) 



For any a / 0, the same factorization method as used 



in Sect. VB can be used, leading to a Casimir force per 
unit surface equal to 



F C /A 



;i + f)r 

8na\L x 

„2r2 7 



[U 3 (e- 2koL *) + 2k L x U 2 (e 



~2k Q L, 



) 



r2k'^LlU 1 (e- 2koL -) - Li 3 (e- 2fcli *) 

-2k x L x ^ 2 {e~ 2k ^) - 2fc?£|Li 1 (e- 2fcl£ -)] , 

(VI. 14) 

where kn = — and k 2 = — + t 2 — ■ In the limit of infinite 
correlation length we have fc — » and fci - 
which we obtain 



.— , from 



F c /A 



1 



[C(3)-Li 3 (e- 



(VI.15) 



8na\L x 

-2k 1 L x U 2 (e- 2klL *) - 2k 2 1 L 2 x Li 1 (e- 2klL *)] . 

The presented result should be compared with , where 
the same system was studied, but a different answer was 
given. At long distances it decays as in the case of white 
noise (VI. 12 1 with a prefactor (1 + |)/a. 



For a — > U, we obtain the quenched limit of the stress 
tensor at the plates 

2 

-r 



(T M (0)) 







dkk 



II ,, : 7T 



Tl x gZ 



(VI. 16) 

This expression is regularized using Eqs. (A.6I and ( |A.3 ) 
with s = 2, resulting in 



F c /A 



r 



dk 



k 1 



„2uL, 



(e 



2uL x 



2ujL x e 
^1) 



2cjL, 



) 



(VI. 17) 

where w = ^/fc 2 + fcg. This integral can be carried out 
to obtain the Casimir force as 

T k L x 



Fc/A 



4it\ 2 k 2 L x e 2kaL * - 1 ' 
In the limit of infinite correlation length we have 

Fc/A ' 



8ir\ 2 n 2 L x ' 



(VI.18) 



(VI.19) 



and in the limit of small correlation length (koL x 3> 1) 
we obtain 



(VI.20) 



As a final case we consider a noise that is rapidly fluc- 
tuating in time but that is homogeneous in space, de- 
scribed by the correlation (V.18). In this 
not diagonal but is given by 

h -QT^ t 1 -(-!)"*] [1 -(-!)"*] , x e r 
"-rim — ZV L - n nO n nO m nO m n. 

'K^n x m x 

(VI.21) 

The stress on the plates is then given by 



(T ra (0)) = 



2r 



oo 



Att 2 E 



[1 _(_!)».][! 



i 2 + 2( 



(VI.22) 



As in the double summation above only odd values of 
n and m are summed, it can be expressed in terms of the 
Elizalde zeta function over odd numbers, defined as 



Zi(a,P,u},8) 



\ " 



1 



^ (a 2 (2n+l) 2 + f3 2 (2m + l) 2 +uj 2 ) s 

(VI.23) 

which can be written in terms of four Elizalde zeta func- 
tions. Using the asymptotic expansion of the Elizalde 
zeta functions given in Eq. (A. 2 1 with p — 2 and s = 1, 



the Casimir force can be expressed as an infinite sum of 
Bessel functions Kq(x) with different values of x. The 
divergent terms, given by the first term in Eq. (A.2), are 



independent of L x , so the final expression is finite and 
given by 



(n,m)/(0,0) 



Fc/A = -"- x Y. E 



K (2V2k L x Vn 2 + m 2 ) 



"5*0 


2y/2k L x 


-2*0 


2y/2k L x 


+ \K 


2V2k L x 



4 1 4 

(VI.24) 

Two limiting cases can be considered to clarify this 
result. First, in the limit of long distances k^L x ^> 1, the 
Casimir force is given by 



(VI.25) 



whereas in the opposite limit of long correlation length 
koL <C 1, the result is 



Fc/A = log (k L) 



(VI.26) 



In the limit of infinite correlation length this result di- 
verges. 
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VII. TWO-FIELD SYSTEM 

In the two systems we have considered so far (reaction- 
diffusion and liquid crystals), the dynamics is described 
by a Hermitian operator and therefore the potential of 
the method described herein is not fully evident. In this 
section, we build a more complex system, described by a 
model with two fields (which could be temperature and 
concentration, for example) coupled in a nonsymmctric 
way. For simplicity and to be concrete we will consider 
that the fields ipi and ijj2 are scalar, subject to Neumann 
boundary conditions, a nd w ith eigenfunctions described 
by the Fourier modes (V.4). In Fourier space, the dy- 
namic equation is 



d_ N, x 
dt \ip2 



a k 

/3k 7k 



(VII. 1) 



The noises are assumed to be white with different corre- 
lation intensities (allowing one of them to be set equal to 
zero later) and no cross-correlation 

(a(r,t)a(r' J 0)=ri<5(r-r')5(t-i / ), 

(a(r,t)6(r',t')} = 0. (VII.2) 

Finally, the stress tensor is assumed to be isotropic, 
depending only on the fields as 



^xx = KltPl + ^2^2 



(VII.3) 



As the dynamic matrix is non-Hermitian, the left and 
right eigenmodes are different, being given by 



/i,k(r) = 
/2,k(r) = 

51, k(r) = 

52, k(r) = 
with eigenvalues 



cos (fc^x) e ik ii r n 
cos(k x x)e ik ^« 

cos(k x x) e ik ii- r " 



^cos(k x x) e ik ii' r ii 



1 



«k-7k 





gk 

1, 



(VII.4) 



A*l,k 
M2,k 



a k , 
7k- 



(VII.5) 



Using these eigenmodes, the different elements needed 
to compute the Casimir pressure are 



^zk.jq 



t ikj'q 



kq 



-r x 



Ok- 7k 



kq 




(VII.6) 



After simple algebra, the stress tensor on the plates is 
obtained as 



: = 2 V Y. 



Ti{ki + K 2j 3 k /(o; k - 7 k ) 2 
2a k 

2TV^/K ~ 7k) 2 
"k + 7k 

(ri^/(a k -7 k ) 2 + r 2 ) K 2 



27k 



(VII.7) 



which for specific models (that is, specific values of a k , 
/3 k , and 7 k ) could be computed and regularized to obtain 
the Casimir force on the plates. 

To show the kind of results that can be obtained we 
consider the simple reaction-diffusion two-field model 
«k = Ai + Dk 2 , /3 k = A12, and 7 k = A 2 + Dk 2 , with 
noise intensities r x = T and T 2 = [35], representing the 
system 



dt 

d-02 
dt 



-AiV>i + DW 2 ^! + f, 

-\2lp2 + -DV 2, 02 - ^12^1- 



(VII. 



Furthermore, we assume that the stress only depends 
on -02, i.e., T xx — nip 2 . Therefore, any eventual Casimir 
force is produced by the fluctuations of the second field 
which are produced by the coupling with the first field. 
The stress on the plates is finally 



(T ra (0)) 



TkX 



12 



2V 



[(Ai + ^fc 2 )(A 2 + Dk 2 ){\i + A 2 + 2Dk 2 



Assuming that Ai 7^ A2, we can apply partial fraction 
decomposition to obtain 



(TWO)) = 



r«A 2 2 



2DV(Xx - A 2 ) 2 



E 



1 



1 



k 2 + k 2 k 2 + k 



k 2 



^3- 



where k\ = \ X /D, k\ = X 2 /D, and k\ = (Ai + A 2 )/2L>. 
Then, we can perform each infinite sum as in the case of 
scalar white noise to obtain the Casimir force as 



0W0)) = 



-T K \ 2 2 



-2fcii x 



^DL X {\ X - A 2 ) 2 



In 



(1 



-2k 2 L x 



(1-e 



-2k 3 L, 



It is interesting to note that, if A x = and/or A 2 = 0, 
this Casimir force diverges. 



CONCLUSIONS 

In this article, we have developed a formalism to study 
Casimir forces in classical systems out of equilibrium 
based on the stochastic dynamical equations of the sys- 
tem under study. The equilibrium case is recovered as 
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a particular limit where the fluctuation-dissipation the- 
orem is valid. 

In particular, we study the interaction which appears 
between intrusions in a medium subject to any kind 
of noise. The study is restricted to additive noise and 
nonquantum systems; quantum cases and multiplicative 
noises are left for future work. The method only relies 
on the stochastic evolution equation of the field in the 
medium, and information about the interaction between 
the medium and the intrusions, as given by the bound- 
ary conditions of the fields at the surface of the bodies. 
No assumptions are made regarding any characteristic of 
the noise, which could be internal of external, thermal or 
induced, white or colored, and even non-Gaussian. 

This formalism reduces to the classical thermal 
Casimir effect when the medium is subjected to an ad- 
ditive Gaussian white noise with autocorrelation ampli- 
tude Q = ksT(C + C + ) and its dynamics is described by 
a Hermitian operator, as shown in (IV. 13 1. 



We have obtained an exact formula for the Casimir 
force felt by a body (IV. 14), which shows how to de- 



rive the Casimir force for any geometrical configuration 
and noise. Equation (IV. 14) can be used to obtain non- 



equilibirum induced self forces over asymmetric bodies, 
as shown in |38j . It can also provide a numerical tool use- 
ful to evaluate Casimir forces for complicated geometries. 
Along this paper, we have used the formula to obtain 
the force between parallel plates in different media (in a 
reaction-diffusion model and in liquid crystals) under the 
influence of different kinds of Gaussian noises, i.e., white 
noise to recover the thermal case already studied in the 
literature, and noises with nonzero spatial or temporal 
correlation lengths, where different forces appears. 

Finally we have shown an example of the evaluation of 
Casimir forces in a system with non-Hermitian evolution 
dynamics, which was an intractable problem until the 
development of the formalism presented herein. 



and admits the asymptotic expansion, valid for all com- 
plex s, 



Z p (s,ai,oj) 



2tt s cj* 



K S _P | 27TW 



£ £ 



i6 Zp-{0} 



2s-p 

2\ t - 




(A.2) 



where K v (z) is the inhomogeneous Bessel function of the 
second kind or Macdonald function. When p = 2, this 
is called the Chowla-Selberg formula. We are mainly 
interested in the case p = 1 with a > 0, for which 



Z 1 (s,a,uj) 



1 



{a 2 n 2 + ui 2 Y 



r(|)r(«-i) x _ 



T(s)a 
Att s 



y(- 

T(s)a ' \aoj 

v ' n— 1 



K s _i ( 2ttw- 



a 



(A.3) 



When studying Casimir forces between plates with 
Dirichlet boundary conditions, the following series needs 
to be computed 



Yi(s,a,o;) = £ 



(a 2 n 2 + uj 2 y 



(A.4) 



It is straightforward to obtain 

Ea 2 n 2 x - 

(a 2 n 2 +u 2 ) s = ^ 



::iZ (a 2 n 2 +u; 2 y 1 



E 



i 



(A.5) 



(a 2 n 2 + lu 2 ) s ' 

and therefore 

Y\(s, a, uS) — Zi(s — 1, a, u>) — ui 2 Zi(s, a,oj). (A. 6) 



APPENDIX: ELIZALDE ZETA FUNCTION 



The computation of the Casimir forces makes use of 
the asymptotic expansion of the Elizalde zeta function, 
which is defined as 136 



Z p (s,ai,u}) = 



(Al) 
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